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We present nanosecond timescale Vlasov-Fokker-Planck-Maxwell modeling of magnetized plasma 
transport and dynamics in a hohlraum with an applied external magnetic field, under conditions 
similar to recent experiments. Self-consistent modeling of the kinetic electron momentum equation 
allows for a complete treatment of the heat flow equation and Ohm’s Law, including Nernst advection 
of magnetic fields. In addition to showing the prevalence of non-local behavior, we demonstrate 
that effects such as anomalous heat flow are induced by inverse bremsstrahlung heating. We show 
magnetic field amplification up to a factor of 3 from Nernst compression into the hohlraum wall. 

The magnetic field is also expelled towards the hohlraum axis due to Nernst advection faster than 
frozen-in-flux would suggest. Non-locality contributes to the heat flow towards the hohlraum axis 
and results in an augmented Nernst advection mechanism that is included self-consistently through 
kinetic modeling. 


Indirect drive inertial confinement fusion (IGF) is ac¬ 
complished through the compression of a fuel capsule due 
to the exposure of an ablator material to radiation. This 
radiation is a result of laser heating of a hollow-cylinder 
made of a high Z material. At the National Ignition Fa¬ 
cility (NIF), a gold hohlraum is irradiated in a precise ar¬ 
rangement of 351 nm laser beam cones. The walls of the 
hohlraum are heated to an electron temperature of ap¬ 
proximately 3-5 keV with the goal of emitting blackbody¬ 
like radiation of 300-500 eV. The resulting x-ray bath is 
intended for uniform compression of the fuel capsule to 
fusion conditions. [1] 

There has been recent interest in the role of applied 
magnetic fields in high-energy-density plasmas [2-4] for 
inertial fusion energy applications. The Magneto-Inertial 
Fusion Electric Discharge System has been developed in 
order to provide steady state magnetic helds for long 
time-scales relative to the experiments. A recent ex¬ 
periment on the Omega Laser Facility with a 7.5 T ex¬ 
ternal axial magnetic field imposed on an Omega-scale 
hohlraum measured a rise in observed temperature along 
the hohlraum axis [5]. Recent modeling showed that hot 
electrons from laser-plasma interactions [6] can be guided 
through the hohlraum, rather than the capsule, using 
such helds [7]. 

From a complete treatment of Ohm’s Law, it has pre¬ 
viously been shown that electron heat transport can ad- 
vect such magnetic helds through the Nernst effect [8-13] 
in addition to well-known MHD processes like “frozen- 
in-how” and resistive diffusion. Dimensionless numbers 
that compare the ratio of the magnitudes of the Nernst 
term in Ohm’s law to that due to bulk plasma how, 
Rn 1 [10], and the Hall term, ^ 1 [12], suggest 


that Nernst convection should be the dominant mecha¬ 
nism for magnetic held transport in a hohlraum. Such 
a hot and semi-collisional environment is, however, also 
rich in non-equilibrium effects that may complicate the 
magnetic held dynamics. 

Laser heating of the plasma results in steep temper¬ 
ature gradients, typically 0{3 keV/50 ptm). The colli- 
sional mean-free-path of a 3 keV electron is 0(10 /im), 
depending on the plasma density. Since Xmip/L < 100, 
non-local effects can be expected to be important [14]. 
The steep temperature gradients that occur due to the 
intense laser heating in a hohlraum have been shown to 
result in non-local heat how [15, 16]. Gareful considera¬ 
tion of the population of electrons with 2vth < v < Avth 
is required as these carry most of the heat. Addition¬ 
ally, inverse-bremsstrahlung heating of a plasma has been 
shown [17, 18] to not only lead to deviations from classical 
transport, prescribed by Braginskii’s transport equations 
[19], but also new transport terms [20]. Both non-local 
transport and laser heating result in modihcations to 
the shape of the distribution function and therefore non¬ 
equilibrium behavior, which mean that classical trans¬ 
port approximations break down. In order to avoid clas¬ 
sical transport approximations, a kinetic approach is nec¬ 
essary. Kinetic modeling allows for the modeling of mag¬ 
netic held dynamics through a self-consistent and gener¬ 
alized Ohm’s Law that can be derived without distribu¬ 
tion function approximations. 

In this Letter, we demonstrate Vlasov-Fokker-Planck- 
Maxwell (VFPM) simulations of a magnetized, full 
hohlraum-scale plasma including ray-tracing of an 
Omega-like laser conhguration over a nanosecond time- 
scale. The hohlraum is considered without an IGF cap- 
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sule but a gas fill throughout. Note that radiation trans¬ 
port and laser-plasma interactions are neglected in these 
calculations, which may change our results if they were 
included. With the use of IMPACTA [21, 22], we stud¬ 
ied the effect of non-equilibrium electron kinetics on ther¬ 
mal energetic and magnetic field dynamics of a Omega- 
scale hohlraum with an externally imposed 7.5 T mag¬ 
netic field. We found that significant proportions of the 
total heat flow are non-local. Additionally, the presence 
of inverse bremsstrahlung heating resulted in anomalous 
heat flow towards the over-dense plasma of the hohlraum 
wall. Therefore, the diffusive heat flow from the laser- 
heated regions is not an adequate description of the ther¬ 
mal energy dynamics. The heat flows from the laser heat¬ 
ing move the externally imposed magnetic field through 
Nernst advection. To examine the effects of Nernst ad- 
vection in relation to the plasma bulk flow, we show mod¬ 
eling without an electron contribution to the transport 
of magnetic field in Ohm’s Law for comparison. 

We find that magnetic field transport due to Nernst 
flow results in significantly faster field cavitation than 
that is possible via frozen-in-flux. Magnetic field cavita¬ 
tion occurs due to heat flow down the density and tem¬ 
perature gradient, which is shown to be non-local. Re¬ 
tention of the distribution function allows for accurate 
modeling of the magnetic field cavitation because the lo¬ 
cal approximation to the Nernst velocity underestimates 
the true convection velocity by a factor of 2. Nernst flow 
into the over-dense region causes magnetic flux pile-up 
at the walls and results in magnetic field amplification 
by a factor of 3. Magnetic flux pile-up does not occur 
with only plasma bulk flow present as there is a negligi¬ 
ble amount of plasma bulk flow toward the wall from the 
laser heated region. 

The Vlasov-Fokker-Planck (VFP) equation for elec¬ 
trons is solved coupled with Ampere’s and Faraday’s 
Laws and a hydrodynamic ion fluid model to describe 
the plasma. The code we use, Impacta [21, 22] uses a 
Cartesian tensor expansion, with the distribution func¬ 
tion expanded as f{t, r, v) = /o -|- fi • v -|- : vv -f ..., 

where is a unit velocity vector. This expansion 

can be truncated in a collisional plasma, as collisions 
tend to smooth out angular variations in the distribu¬ 
tion function, resulting in a close to isotropic distribu¬ 
tion, represented by /q. Higher orders are successively 
smaller perturbations, ^ fi ^ /o etc. In the classical 
limit that /o is a Maxwell-Boltzmann velocity distribu¬ 
tion, Impacta has been shown to agree with Braginskii’s 
transport equations [21]. In Impacta, terms up to and 
including can be retained. These simulations, however, 
are collisional enough such that may be neglected to 
an error 0(Anifp/L)^. 

A 2-dimensional slice of a hohlraum is modeled in 
the x-y plane where the y-axis represents the longitu¬ 
dinal axis of the hohlraum and the fuel pellet would 
sit at the origin. The hohlraum walls are represented 
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FIG. 1: (a) Ray tracing profile overlaid onto laser 
intensity profile (W/cm^) at t = 0. (b) Electron Plasma 

Temperature (keV), (c) Externally applied magnetic 
field (T), (d) Nernst Velocity (uAr/r'tho) at t = 250 ps. 

by a dense, high-Z plasma located at approximately 
X = ±800 /rm, and the gas fill by a low-Z plasma, 
with the overall Z distribution described by the function 
Z{x,y) = 59.25 ± 19.75 tanh(^^^). Electron number 
density is described by the function n^{x,y) = (2.98 ± 
2.93tanh(2^)) x lO^^ cm-^. The initial uniform tem¬ 
perature was ksTeo = 160 eV. The initial uniform mag¬ 
netic field was Bo(y) = 7.5 T and InAg^ = 5.4. To con¬ 
vert from the normalized units, rieo = 5 x 10^° cm“^ 
and Utho/c = 0.025 are used. . The laser parameters are 
designed to resemble those of ref. [5]. The ray tracing 
package tracks the three beam cones that enter at 21, 42, 
and 59 deg from the axis, to their respective refraction 
points and allows for some reflection. 

The rays and the initial heating profile are shown in 
fig. la. Eigure lb shows the temperature profile after 
300 ps of laser heating. Eigure Ic shows the cavitation 
and amplification in the in-plane magnetic field profile 
caused by intense laser heating. The Nernst velocity, 
shown in fig. Id, is directed towards the hohlraum axis 
in the low density gas hll and into the hohlraum wall in 
the Au plasma. Throughout the rest of this Letter, we 
show that the Nernst flow is primarily responsible for the 
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magnetic field profile seen in fig. Ic. 

It should also be noted that while self-generated mag¬ 
netic fields, Bz in this geometry, are also included in this 
study, their dynamics are not the subject of this Letter. 
Many of the effects discussed here, however, can be ap¬ 
plied to the self-generated fields and will be described in 
detail separately. 
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FIG. 2: (a) Heat flow (weneot't^o) (b) 1 - 9eq. (i)/'?code 
(c) V^{fcode - /mb) (d) U®(/code “ /sg) 

@ X = 0.4 mm, y = —0.6 mm, 
m = 2.625 , t = 100 ps. 

Inverse bremsstrahlung heating of the plasma results 
in a super-Gaussian electron distribution [17], which con¬ 
sequently modifies the transport coefficients [20, 23] and 
even introduces new terms including an anomalous heat 
flux up a density gradient q„, represented by the last 
term in eq. (1); 

Qe = ■ j “ (^ + ’ (1) 

where ^/>, (p and n are transport coefficients as described 
in reference [20]. q„ increases as m > 2 increases, where 
m is the power of the super-Gaussian distribution func¬ 
tion defined by/ sg{?^) = C'(m)ne/uthexp(- (v/aeVth)"^) 
where ae = [3r(3/m)/2r(5/m)]^/^ and C{m) = 

m/4TTalT{3/m). 

In these simulations, by finding the best fit of a super- 
Gaussian distribution to /o, to reaches a maximum of 3.1 
near the centers of the laser heated regions, but varies 
spatially and temporally, thus requiring the preservation 


of the distribution function at each point throughout the 
simulation for accurate calculation of the heat flow. Us¬ 
ing the theory detailed in refs. [20, 23], the heat flow 
can be modified in a hydrodynamics code to include this 
effect. However, the distribution is not precisely a super- 
Gaussian due to other effects such as non-locality and 
therefore this fix remains an approximation. 

We examine the relative magnitudes of the real heat 
flow, and the classical heat flow calculated using all three 
terms that form the full post-processed heat flow from 
eq. (1) that includes anomalous heat flow. Calculation 
of the anomalous heat flow as a function of the best-fit 
distribution function, table look-up, and pressure gradi¬ 
ent shows that there is heat flow into the hohlraum wall 
due to the pVPf. term and this approximately results in 
a 10% correction to the diffusive heat flow i.e. kVT^. 

A majority of the disagreement between the heat flow 
from the code and the heat flow from the post-processed 
modified classical transport theory is due to the strongly 
non-local heat flow that is prevalent in the hohlraum. 
Figure 2b shows a 2D profile of a metric for quantifying 
the magnitude of the discrepancy between the two heat 
flows, described by the relative difference between the 
classical and calculated heat fluxes, 1 — geq. (i)/'Zcode- 

The regions within the black contours have ±25% 
agreement between the two heat flows. The white con¬ 
tours correspond to regions of high non-locality where 
the classical transport calculation is an underapproxima¬ 
tion, while the blue contours correspond to regions where 
the heat flow is significantly overcalculated by classical 
transport. Heat flow from regions near the temperature 
hotspots, ±50 /rm, is overestimated by the classical calcu¬ 
lation while the heat flow further away from the hotspots, 
±200 /TTO, is underestimated, as expected from the exis¬ 
tence of non-locality. The regions of relative agreement 
are ± 50 — 200/rm from the hot spots. Due to the laser 
heating, the thermal electron mean-free-path increases, 
Amfp/L > 0.02, suggesting that non-local heat flow be¬ 
comes prevalent in the laser heated region. 

Consideration of the in-plane electron distribution 
function f{9,v) = /o ± fixVx ± fiyVy can show the 
significance of inverse-bremsstrahlung heating and non¬ 
locality. Since q oc J f {6, v)^{9, (j))dv sin 9d6dp, the 
important contributions to the heat flow may be best 
illustrated by the function v^f(9,v). Figure 2c and d 
show the difference between the calculated distribution 
v^f and (c) a Maxwell-Boltzmann u®/mb and (d) a super- 
Gaussian with best fit to to, both with Tg equal to that 
of f{x = 0.4,?/ = —0.6). Figure 2c shows that / > /mb 
in the region 2 < uth < 4 and / < /mb in the re¬ 
gion 4 < vth < 6, which is characteristic of inverse- 
bremsstrahlung heating. Calculating the heat flow con¬ 
tribution difference between the real distribution and the 
best-fit super-Gaussian (to « 2.625 in this case), shows 
that the inverse-bremsstrahlung model does not replicate 
the distribution function fully due to anisotropy from 
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the flow and non-local effects. The enhanced tail and 
shifted center in the 180° direction is characteristic of 
the (non-local) heat flow down the density gradient while 
the colder return flow is a result of the features in the 0° 
direction. 

As shown in ref. [9], the Nernst velocity is, 


Vat 


M J_ 
2(u3) erif. 
K ■ VTe 

5/2Pe ’ 


( 2 ) 

(3) 


It can be shown for this geometry that By has no field 
generation terms from the curl of Ohm’s Law and there¬ 
fore, can be transported through (v^r -I- C) x B term in 
addition to resistive diffusion. Over 0.5 ns, the simulation 
shows that there is magnetic field cavitation resulting in 
flux pile-up on the hohlraum axis and significant com¬ 
pression at the hohlraum wall due to the energy deposi¬ 
tion from the laser. Pile-up of the magnetic flux results 
in a 25 T magnetic field, more than 3 times the strength 
of the initial 7.5 T field. 



FIG. 3: Magnetic field (T) after 50 ps with only plasma 
bulk flow (a) and full Ohm’s Law (b). 

Magnetic field after 400 ps with only plasma bulk flow 
(c) and full Ohm’s Law (d) 


In order to determine the effect of Nernst advection 


on the magnetic field evolution, simulations with and 
without the B x fi term in the fi equation were com¬ 
pared. This term is responsible for the interaction of 
kinetic electrons with the magnetic field. It is respon¬ 
sible for the Nernst and Hall terms in Ohm’s Law as 
well as the Righi-Leduc effect in the heat flow equation. 
Simulations agree with the previous determination that 
j vtv because <C 1 and therefore, the Hall effect 
can be neglected. The magnetic field after 50 ps without 
and with full Ohm’s Law treatment is shown in fig. 3a 
and fig. 3b, respectively. The laser heated region results 
in magnetic field cavitation in both cases but the magni¬ 
tudes differ. It is not evident in fig. 3a since the field is 
only modified by a few percent by the plasma bulk flow. 
Thermal energy transport results in a more noticeable 
change immediately over 50 ps. 

An estimate of the time-scale for the plasma bulk flow 
to transport frozen-in magnetic fields to the center of the 
hohlraum is given by, ^ « -j=M== ~ 2 ns. Figure 3d 

shows that including the Nernst effect results in mag¬ 
netic field cavitation on a faster time-scale than can be 
expected due to field advection only through bulk plasma 
flow in fig. 3c. In the case of a 7.5 T initial field strength, 
the magnetic field on the axis grows to 30 T within 0.5 
ns. Figure 3d also shows that the magnetic flux pile-up in 
the hohlraum wall occurs due to the Nernst effect. The 
field increases to a strength of nearly 25 T towards the 
hohlraum wall. 



FIG. 4: The magnetic field at the hohlraum axis 
decreases as applied field strength is increased because 
the Nernst effect is mitigated at higher magnetic field 
strengths, preventing the magnetic field from 
accumulating near the hohlraum axis, (t = 300 ps). 
The discrepancy between the approximated and exact 
Nernst velocity also decreases. 


We also ran a series of simulations with varying initial 
applied field ByQ to understand how the field strength af¬ 
fects the hohlraum dynamics. Figure 4 shows the results 
of the magnetic field cavitation study for increasing field 
strengths suggesting that as 

lim By—o^^is/ByQ = 1. (4) 

ByQ—>-(X> 

The maximum value of in the domain of magnetic 
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field advection towards the axis (—0.5 mm < x < 
0.5 mm) is chosen. This trend can be explained by the 
observed reduction in the Nernst velocity towards the 
hohlraum axis as the magnetization increases (also shown 
in fig. 4), which quenches magnetic field transport. These 
V]\[{ujt) curves are in line with other predictions [9, 12, 20] 
that vn oc l/wr for uir ^ 1. Figure 4 also shows that 
the exact Nernst velocity from eq. (2) is consistently, and 
up to 2 X larger than what the local approximation from 
eq. (3) would predict for T^, rig, and B profiles at 300 ps. 
This discrepancy decreases at higher field strengths due 
to magnetic field induced localization of the heat flow 
carrying electrons. 

The degree of magnetic flux pile up in the hohlraum 
wall, however, is not so strongly affected by the increase 
in magnetic field strength because lot ~ n“^. The mag¬ 
nitude of maximum field strength in the wall ranges from 
2 < By/Byo < 3 for 1 < Byo < 100 T. 

We have shown Vlasov-Fokker-Planck modeling of an 
external magnetic field of 1-100 T imposed upon a 
Omega-scale hohlraum. Magnetic flux pile-up causes an 
increase in magnetic field magnitude by a factor of 3 
for a 7.5 T magnetic field. Additionally, the heat flow 
is responsible for magnetic field cavitation on a faster 
time-scale than that from the bulk flow of the plasma. 
Not only is the heat flow strongly non-local, it also 
has distinct signatures of inverse bremsstrahlung heat¬ 
ing. The ability to preserve distribution function infor¬ 
mation through use of a kinetic code allows to model the 
heat flow accurately. Full Vlasov-Fokker-Planck-Maxwell 
treatment of the system enables accurate modeling of 
magnetic field dynamics. We have shown that the Nernst 
flow is the dominant mechanism for magnetic field trans¬ 
port and is responsible for the increase in field strength, 
up to 100 T for a initial 100 T field, in the wall as well 
as cavitation of the magnetic field towards the hohlraum 
axis. The magnetic field cavitation is mitigated at higher 
field strengths. Furthermore, the Nernst velocity is up to 
2x larger in VFPM than would be predicted by classical 
transport. 

These findings suggest that attempting the same cal¬ 
culation with a classical description of transport would 
result in significantly different B & Tg evolution. Accu¬ 
rate modeling of these quantities has implications for con¬ 
trolling levels of laser plasma interactions [5, 6] and hot 
electron propagation [7] in the gas fill and understand¬ 
ing the hot spots on the dense wall that generate X-rays. 
The enhanced electron transport and B field physics pre¬ 
sented here could affect details of X-ray drive if incorpo¬ 
rated into full-scale radiation-hydrodynamics modeling 
(including reduced phenomenological laser-plasma inter¬ 
action models) of indirect drive with externally applied 
B-field. 

The authors would like to thank A. Hazi, J. Moody, D. 
Strozzi of LLNL and A. Sefkow of Sandia for useful dis¬ 
cussions regarding related work at NIF and Omega. The 


modeling was performed using computational resources 
and services provided by Advanced Research Comput¬ 
ing at the University of Michigan, Ann Arbor. This re¬ 
search was supported by the DOE through Grant No. 
DE SC0010621. 


* archisj@umich.edu 
^ agrt@umich.edu 

[1] J. D. Lindl, P. Amendt, R. L. Berger, S. G. Glendinning, 
S. H. Glenzer, S. W. Haan, R. L. Kauffman, O. L. Lan- 
den, and L. J. Suter, Physics of Plasmas 11, 339 (2004). 

[2] P. Y. Chang, G. Fiksel, M. Hohenberger, J. P. Knauer, 
R. Betti, F. J. Marshall, D. D. Meyerhofer, F. H. Seguin, 
and R. D. Petrasso, Physical Review Letters 107 , 035006 
( 2011 ). 

[3] M. Hohenberger, P.-Y. Chang, G. Fiksel, J. P. Knauer, 

R. Betti, F. J. Marshall, D. D. Meyerhofer, F. H. Seguin, 
and R. D. Petrasso, Physics of Plasmas 19 , 056306 
( 2012 ). 

[4] G. Fiksel, W. Fox, a. Bhattacharjee, D. Barnak, P.-Y. 
Chang, K. Germaschewski, S. Hu, and P. Nilson, Phys¬ 
ical Review Letters 113 , 1 (2014). 

[5] D. S. Montgomery, B. J. Albright, D. H. Barnak, P. Y. 
Chang, J. R. Davies, G. Fiksel, D. H. Froula, J. L. Kline, 
M. J. MacDonald, a. B. Sefkow, L. Yin, and R. Betti, 
Physics of Plasmas 22 , 010703 (2015). 

[6] S. P. Regan, N. B. Meezan, L. J. Suter, D. J. Strozzi, 
W. L. Kruer, D. Meeker, S. H. Glenzer, W. Seka, 

C. Stoeckl, V. Y. Glebov, T. C. Sangster, D. D. Meyer¬ 
hofer, R. L. McCrory, E. a. Williams, O. S. Jones, D. a. 
Callahan, M. D. Rosen, O. L. Landen, C. Sorce, and 
B. J. MacGowan, Physics of Plasmas 17 , 020703 (2010). 

[7] D. J. Strozzi, L. J. Perkins, M. M. Marinak, D. J. 
Larson, J. M. Koning, and B. G. Logan, (2015), 
arXiv: 1508.00803. 

[8] A. Nishiguchi, T. Yabe, M. G. Haines, M. Psimopou- 
los, and H. Takewaki, Physical Review Letters 53, 262 
(1984). 

[9] M. G. Haines, Plasma Physics and Controlled Fusion 28 , 
1705 (1986). 

[10] C. P. Ridgers, R. J. Kingham, and A. G. R. Thomas, 
Physical Review Letters 100 , 075003 (2008). 

[11] L. Willingale, a. G. R. Thomas, P. M. Nilson, M. C. 
Kaluza, S. Bandyopadhyay, a. E. Dangor, R. G. Evans, 
P. Fernandes, M. G. Haines, C. Kamperidis, R. J. King- 
ham, S. Minardi, M. Notley, C. P. Ridgers, W. Rozmus, 
M. Sherlock, M. Tatarakis, M. S. Wei, Z. Najmudin, and 
K. Krushelnick, Physical Review Letters 105, 1 (2010). 

[12] A. S. Joglekar, A. G. R. Thomas, W. Fox, and A. Bhat¬ 
tacharjee, Physical Review Letters 112 , 105004 (2014). 

[13] C. Li, F. Seguin, J. Frenje, N. Sinenian, M. Rosen¬ 
berg, M.-E. Manuel, H. Rinderknecht, A. Zylstra, R. Pe¬ 
trasso, P. Amendt, O. Landen, A. Mackinnon, R. Town, 

S. Wilks, R. Betti, D. Meyerhofer, J. Soures, J. Hund, 
J. Kilkenny, and A. Nikroo, Nuclear Fusion 53 , 073022 
(2013). 

[14] D. R. Gray, J. D. Kilkenny, M. S. White, P. Blyth, and 

D. Hull, Physical Review Letters 39, 1270 (1977). 

[15] G. Gregori, S. H. Glenzer, J. Knight, C. Niemann, 
D. Price, D. H. Froula, M. J. Edwards, R. P. J. Town, 



6 


A. Brantov, W. Rozmus, and V. Y. Bychenkov, Physical 
Review Letters 92 , 205006 (2004). 

[16] J. Hawreliak, D. M. Chambers, S. H. Glenzer, A. Gou- 
veia, R. J. Kingham, R. S. Marjoribanks, P. a. Pinto, 
O. Renner, P. Soundhauss, S. Topping, E. Wolfrum, P. E. 
Young, and J. S. Wark, Journal of Physics B: Atomic, 
Molecular and Optical Physics 37 , 1541 (2004). 

[17] A. B. Langdon, Physical Review Letters 44, 575 (1980). 

[18] J. M. Liu, J. S. De Groot, J. P. Matte, T. W. John¬ 
ston, and R. P. Drake, Physical Review Letters 72 , 2717 


(1994). 

[19] S. 1. Braginskii, Reviews of Plasma Physics 1, 205 (1965). 

[20] C. P. Ridgers, a. G. R. Thomas, R. J. Kingham, and 
a. P. L. Robinson, Physics of Plasmas 15, 092311 (2008). 

[21] R. Kingham and A. Bell, Journal of Gomputational 
Physics 194 , 1 (2004). 

[22] a. G. R. Thomas, R. J. Kingham, and G. P. Ridgers, 
New Journal of Physics 11 , 033001 (2009). 

[23] J. J. Bissell, C. P. Ridgers, and R. J. Kingham, New 
Journal of Physics 15 , 025017 (2013). 



